Journal  of  Power  Sot 


;  212  (2012)  226-232 


Contents  lists  available  at  SciVerse  ScienceDirect 

Journal  of  Power  Sources 


ELSEVIER 


epage:  www.e 


■n/locate/jpowsou 


Short  communication 


Sliding  mode  control  of  a  linearized  polymer  electrolyte  membrane  fuel 
cell  model 


Gunhyung  Park*,  Zoran  Gajic 

Department  of  Electrical  and  Computer  Engineering,  Rutgers  University,  94  Brett  Road,  Piscataway,  NJ  08854-8058,  USA 


ARTICLE 


N  F  O 


A  B  S  T  R 


C  T 


Article  history: 

Received  13  December  2011 
Received  in  revised  form 
3  April  2012 
Accepted  5  April  2012 
Available  online  13  April  2012 


Keywords: 

PEM  fuel  cell 
Model  linearization 
Sliding  mode  control 


The  well-known  nonlinear  fifth-order  model  of  a  proton  exchange  membrane  (PEM,  also  known  as 
polymer  electrolyte  membrane)  fuel  cell  (PEMFC)  appears  to  be  pretty  complex.  In  this  paper,  we 
introduce  the  linearized  model  of  the  original  nonlinear  system  and  propose  a  sliding  mode  technique  to 
keep  the  pressures  of  hydrogen  and  oxygen  at  the  desired  values  despite  of  changes  of  the  fuel  cell 
current.  Since  the  equilibrium  point  at  steady  state  is  unique,  we  perform  Jacobian  linearization  of  the 
original  model  at  steady  state  and  find  the  state  space  matrices  of  the  linearized  system  using  the 
MATLAB  Symbolic  Tool  Box.  The  linearized  system  is  asymptotically  stable  as  well  as  controllable  and 
observable.  In  this  paper  we  show  that  a  sliding  mode  control  technique  copes  very  well  with  the  fuel 
cell  external  disturbance  (changes  of  the  fuel  cell  current)  and  produces  excellent  results  for  hydrogen 
and  oxygen  required  pressures. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  electrochemical  energy  devices  that  convert  the 
chemical  energy,  during  a  hydrogen-oxygen  reaction,  into  elec¬ 
tricity,  heat,  and  water.  As  a  renewable  energy  source,  fuel  cells  are 
one  of  the  promising  energy  technologies  with  high  efficiency  and 
low  environmental  impact.  Proton  exchange  membrane  fuel  cells 
are  the  most  developed  and  popular  type  of  fuel  cells,  using 
hydrogen  as  the  fuel.  PEMFC  represents  a  nonlinear,  multiple-input 
and  multiple-output  dynamic  system  [1], 

Third-order  models  of  PEMFC  can  be  found  in  [2]  (linear  model) 
[3]  (bilinear  model),  and  [4]  (nonlinear  model).  Na  and  Gou  have 
derived  a  fifth-order  nonlinear  model  [5]  since  Chiu  et  al.  [4] 
indicated  a  need  for  using  higher-order  dimensional  models  of 
anode  water  (needed  for  membrane  humidification)  and  cathode 
nitrogen,  see  Fig.  1.  A  nonlinear  ninth-order  model  of  PEMFC  was 
derived  in  [6],  The  model  of  [6]  and  its  simplified  fourth  and  sixth- 
order  variants  were  considered  in  [7,8], 

In  the  last  few  years,  several  control  strategies  for  PEMFC  are 
proposed.  The  sliding  mode  control  technique  that  is  robust  against 
disturbances  has  been  recently  considered  in  several  papers  [7-15], 
A  second-order  sliding  mode  controller  is  designed  for  the 
breathing  subsystem  of  a  PEMFC  stack  in  [8],  where  the  authors 
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have  focused  on  elimination  of  the  chattering  phenomenon.  The 
model  used  in  [8]  is  a  nonlinear  model  of  order  six,  which  is  derived 
based  on  the  work  of  Pukrushpan  et  al.  [6],  Talj  et  al.  [7],  have  first 
simplified  and  reduced  the  ninth-order  model  of  [6]  to  a  fourth- 
order  highly  nonlinear  model  and  experimentally  justified  such 
a  procedure.  Then,  they  designed  the  corresponding  sliding  mode 
controller  using  as  the  sliding  variable  the  difference  between  the 
actual  and  nominal  angular  air  compressor  speeds.  The  oxygen  flow 
problem  with  real  time  implementation  of  a  sliding  mode  controller 
has  been  considered  for  the  first-order  model  that  is  obtained  from 
the  process  input/output  data  in  [12],  In  [13],  a  hybrid  controller 
composed  of  an  internal  mode  control  based  PID  controller  and  an 
adaptive  sliding  mode  controller  has  been  designed.  The  first 
controller  is  used  to  control  the  hydrogen  reformer  and  the  second 
controller  is  used  to  control  the  PEMFC  model  based  on  the  work  of 
[2],  A  sliding  mode  control  scheme  is  proposed  for  the  DC/DC  buck 
converter  that  guarantees  a  low  and  stable  output  voltage  given 
transient  variations  in  the  output  voltage  of  a  PEMFC  in  [  14],  A  fuzzy 
sliding  mode  current  controller  of  a  hybrid  fuel  cell/energy-storage 
systems  is  considered  in  [10].  The  method  is  presented  for  designing 
controllers  for  DC/DC  and  DC/AC  converters.  In  [11],  Hajizadeh  and 
his  coworkers  used  the  fuel  cell  model  of  [11  ]  coupled  with  a  simple 
second-order  model  for  the  hydrogen  reformer  and  a  linear  super 
capacitor  model  to  design  a  sliding  mode  controller  for  active  power 
under  unbalanced  voltage  sag  conditions.  A  sliding  mode  controller 
of  DC/DC  converters  for  a  simplified  dynamic  model  for  fuel  cells  is 
used  in  [9],  Li  et  al.  have  presented  a  rapid-convergent  sliding  mode 
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Fig.  1.  Na  and  Gou  2008  model  Schematic  [5], 

controller  for  the  temperature  control  system  of  PEMFC  stack  [15]. 
That  paper  did  not  consider  the  fuel  cell  connection  to  an  electric 
grid  [9—11,14],  but  the  model  has  included  state  variables  that 
represent  the  fuel  cell  temperature  [15]. 

In  this  paper  we  propose  a  sliding  mode  controller  design  for 
the  fifth-order  nonlinear  model  of  PEMFC  developed  in  [5],  see 
also  Ref.  [16].  The  state  variables  in  this  model  represent  respec¬ 
tively  the  pressures  of  hydrogen  and  water  at  the  anode  side  and 
the  pressures  of  oxygen,  nitrogen,  and  water  at  the  cathode  side, 
that  is 

X(t)  =  [PH2(t)  PH2oJt)  Po2(t)  PNz(t)  PH20c(t)]T  (!) 

=  [xi(t)  x2(t)  x3(t)  x4(t)  x5(t)]T 
The  state  space  model  is  given  by 
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where  R  is  the  universal  gas  constant,  T  is  temperature,  V4,  \ fc  and 
are  anode  and  cathode  volumes.  Ci,  C2  are  known  constants  [5,16], 
cpa,  <pc  are  the  relative  humidity  constants,  P„s  is  the 
saturation  pressure,  /.H2,  /air  and  are  stoichiometric  constants, 
Yh2  =  0.99, y02  =  0.21, Yn2  =  0.79,  are  reactant  fractions.  I  is  the 
cell  current,  and  it  is  considered  as  a  disturbance  since  it  changes  as 
Vfc/Ri,  where  VfC  is  the  produced  fuel  cell  voltage  and  Rl  is  the  load  of 
active  users,  which  changes  randomly,  y  is  the  coefficient  measured 
experimentally  in  [17]  for  the  back-diffusion  of  water  from  the 

cathode  to  the  anode  which  is  defined  as  H2Oback  =  YH2Ombr  where 

ThOback  is  the  water  back-diffusion  flow  rate  and  H2Ombr  is  the 
membrane  water  inlet  flow  rate  [  16],  In  this  model,  it  is  assumed  that 
the  water  in  the  system  is  perfectly  controlled  and  the  back-diffusion 
is  not  considered  which  means  y  is  zero. 

The  model  output  variables  are  defined  by  [5,16] 


The  system  control  input  is  given  by  u(t)  —  [ua(t)  uc(t)]T 
where 


ua(t)  =  l(H2in(t)  +  H20Ain(t))  (4) 

H2m  (t)  and  H2OAin  represent  inlet  flow  rates  of  the  anode  side 
hydrogen  and  water  vapor  with  ka  being  a  known  constant,  and 

Uc(D  =  jg  <02in(t)  +  N ato(t)  +  H2Ocin(t))  (5) 

02in  ( t ),  N2in  (t)  and  H2Oan  (t)  represent  respectively  inlet  flow 
rates  of  the  cathode  oxygen,  nitrogen,  and  water  with  kc  being 
a  known  constant. 


2.  Linearization  of  PEM  fuel  cell  dynamic  model 


Using  MATLAB  Symbolic  Math  Tool  box,  we  obtain  a  unique 
equilibrium  point  for  the  given  system  (2)  as 
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Functions  ipy  are  given  in  Appendix  A. 

Using  the  Jacobian  linearization  technique  [18,19],  the  system 
(2)  can  be  linearized  at  the  equilibrium  point.  The  Jacobian  matrices 
at  the  equilibrium  point  defined  by  x,  u,  and  I,  corresponding  to 
system  (2),  and  represented  in  general  as  x  =  f(x,u,I),  are 


m  m  fj 

9s|x=Xju=BjJ=j’  0u|x=Xju=u,/=/’  07lx=x,u=u,/=l 

where 


(7) 


x(t)  =x+  dx(t) 

u(t)  =  u  +  6u(t)  (8) 

x(I)  =  7 +  5,(1) 

The  perturbations  defined  in  (8)  are  assumed  to  be  small  [18,19], 
The  linearized  system  is  defined  by 


6 x(t)  =A6x(t)  +  Bdu(t)  +  G6I(t) 
with  the  constant  matrices  given  by 
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It  has  been  found  that  the  matrix  A  is  given  by 
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Elements  a, j  are  given  in  Appendix  B. 


(12) 
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The  matrix  B  is  similarly  obtained  a: 


b32 
b42 
b52 

with  elements  by  given  in  Appendix  C. 
The  matrix  G  is  obtained  as  follows 


3.  Sliding  mode  controller  d 
cell  dynamic  model 


(14) 

n  of  the  linearized  PEM  fuel 


Sliding  mode  control  is  a  form  of  variable  structure  control  [20], 
where  sliding  surfaces  are  designed  such  that  systems  trajectories 
exhibit  desirable  properties.  A  system  using  sliding  mode  control 
has  been  considered  as  a  robust  system,  which  yields  to  reduced 
system  sensitivity  to  uncertainties  and  exogenous  disturbances. 

Sliding  mode  control  systems  have  been  studied  in  different  set¬ 
ups  by  many  researchers  [20],  The  controller  is  designed  in  two 
steps  —  finding  the  sliding  surface  and  reaching  the  sliding  mode. 
After  finding  sliding  surfaces  using  the  method  of  [20]  for  linear 
systems  or  the  Lyapunov  method  for  nonlinear  systems  [21],  the 
design  of  sliding  mode  control  is  achieved  as  follows.  Firstly,  a  sliding 
surface  is  defined  which  ensures  that  the  system  remains  on  a  hyper¬ 
plane  after  reaching  it  from  any  initial  condition  in  a  finite  time. 
Secondly,  discontinuous  control  is  designed  to  render  a  sliding  mode. 
Approaches  [20-23]  can  be  used  for  continuous-time  sliding  mode 
control  which  has  been  recognized  as  a  robust  control  approach, 
which  yields  to  reject  matched  disturbances  and  system  uncer¬ 
tainties.  The  matching  condition  [22],  provided  the  control  input 
makes  the  system  asymptotically  stable,  assures  robustness  against 
parametric  uncertainties  and  exogenous  disturbances. 

In  the  following,  we  present  some  basic  concepts  of  the  sliding 
mode  control  technique  for  linear  systems  such  as  the  sliding 
surface  design,  sliding  mode  control  design,  and  the  disturbance 
rejection  condition  [23]. 

Consider  a  continuous-time  linear  system 

x(t)  =  Ax(t)  +  Bu(t)  (15) 

where  x(t)  e  Rn,  u(t)  e  Rm,  and  A,  B  are  constant  matrices  of 
appropriate  dimensions,  and  B  has  full  rank.  There  exists  a  simi¬ 
larity  transformation  defined  by  [20] 


and  columns  of  the  n  x  (n  -  m)  matrix  N  composed  of  basis  vectors 
in  the  null  space  of  Br,  which  puts  Equation  (15)  into  the  form 


q(t)  ^  Aq(t)  +  Bu(t)  (18) 

with  A  =  HAH  1  andB  =  HB  =  [g  ].  Equation  (18)  is  decomposed 
as  follows 


r  9i(t)i  _  r^n  A12]  . 

[  92(f)  J  1^21  A22J  [q2(t)\ 


£ u(t) 


(19) 


where  (ji(t)eRn_m,  q2 (t)eRm,  and  Br  is  an  m  x  m  nonsingular 
matrix. 

Equation  (19)  yields 


9i(t)  =  An<?i(t)  +Ai2<?2(t)  (20) 

and 

92(f)  =  A2i9i(f)  +A22q2(t) +Bru(t)  (21) 

q2 (t)  is  treated  as  a  control  input  to  the  system  (20)  and  a  state 
feedback  gain  K,  which  makes  the  system  asymptotically  stable,  is 
defined  by 


92(f)  =  -Kq,(t).  (22) 

For  the  system  (20),  Utkin  and  Young  have  shown  that  (An ,  A12)  is 
controllable  if  and  only  if  (A,  B )  is  controllable  [20],  see  also  Ref.  [24], 
On  the  sliding  surface,  the  system  trajectory  in  the  (q  1  (t) ,  q2  (t) ) 
coordinates  is  expressed  as 


I*  «[£$]  -°  <23> 

or 

s(t)  =  [/<■  Im]Hx(t)  =  Gx(t)  =  0  (24) 

in  the  original  coordinates.  Consider  the  same  system  as  in  (15) 
with  a  disturbance  d(t) 

x(t)  =  Ax(t)  +  Bu(t)+Ed(t)  (25) 

The  sliding  variable  dynamics  is  given  by  (the  sliding  surface  is 
chosen  as  s(t)  =  Gx(t)) 

s(t)  =  Gx((§  =  CAx(t)  +  GBu(t)  +  GEd(t)  (26) 

The  control  law  which  satisfies  the  reaching  condition  directly 
can  be  chosen  as 

u(t)  =  -(CB^CMt)  ~  (GB)-1  CT  +  «)  W 

where 

r  =  l|CE||dmax  (28) 

The  disturbance  matching  condition  [22],  is  given  by 


9(f)  =  Hx(t) 

(16) 

with 

H  =  [N  B]r 

(17) 

rank([B  £])  =  rank([B])  (29) 

For  the  state  space  model  (2),  our  objective  is  to  keep  yref  in 
a  certain  range,  which  means  to  keep  e  =  y  —  yref  around  0.  We  can 
define  a  sliding  surface  as  follows 
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(30)  The  corresponding  linearized  system  is 


which  yields 


s  =  y  -  yref  =  y  = 


(31) 


The  sliding  mode  control  laws  that  satisfy  ss  <  0  can  be  deter¬ 
mined  from  Equations  (32)  and  (33),  will  be  presented  in  the  next 
section. 

Several  other  more  complicated  techniques  for  designing  sliding 
surfaces  and  sliding  mode  controllers  can  be  found  in  the  engi¬ 
neering  literature,  see  for  example  Refs.  [25-27],  However,  they  are 
not  needed  for  the  purpose  of  this  paper,  since  with  the  already 
introduced  two  standard  and  simple  sliding  mode  techniques  we 
get  excellent  results.  As  it  will  be  shown  in  the  simulation  results, 
the  chattering  phenomenon  (the  main  problem  in  using  sliding 
mode  control)  will  be  fully  suppressed  and  the  abrupt  changes  of 
the  fuel  cell  current  (disturbance  in  a  pretty  broad  range  of 
80  A— 200  A)  will  have  no  impact  on  hydrogen  and  oxygen 
pressures. 


4.  Linearized  PEM  fuel  cell  dynamic  model 
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By  examining  the  eigenvalues  of  the  matrix  A  we  found  that  all  of 
them  are  in  the  left  half  complex  plane,  so  that  this  system  is  open- 
loop  asymptotically  stable.  The  output  equation  according  to  Equa¬ 
tion  (3)  has  been  already  defined  in  the  linear  form  in  [5]  and  [16]  as 


4.1.  Model  linearization 

The  numerical  data  taken  from  [5],  used  in  this  paper,  are  pre¬ 
sented  in  Table  1. 

With  the  realistic  numerical  data  in  Table  1,  we  have  found  that 
the  considered  fuel  cell  system  has  a  unique  equilibrium  point 
(steady  state  point)  given  by 


y(t)  =  Cx(t) 

with 


[1  0  0  0  01 

0  0  10  0 


(34) 


(35) 


x(t)  =  [2.6509  0.0003  7.009  26.175  0.3390]T  (32) 


Having  obtained  matrices  A,  B,  and  C,  we  can  examine 
the  controllability  and  observability  of  the  linearized  system  [28], 
The  controllability  can  be  examined  by  studying  the  rank  of  the 
controllability  matrix  defined  by  [28] 


Tablet  Co  =  [b  AB  A2B  A3B  A4b\  (36) 

Parameters  of  fuel  cell.  L  J 


Symbol  Parameter 


Va 

kc 


H21n 

H20Ain 

02i„ 

N2i„ 

H2Oci„ 

I 


Gas  constant 

Operating  cell  temperature 
Number  of  cells 

Cathode  volume 
Anode  conversion  factor 
Cathode  conversion  factor 
Fuel  cell  active  area 
Faraday  constant 
Saturation  pressure 
02  reactant  factor 

H2  reactant  factor 
N-A/2F 
1.2684 NA/F 

H2  stoichiometric  constant 
Air  stoichiometric  constant 
Anode  humidity  constant 
Cathode  humidity  constant 
H2  inlet  flow  rate 
H20  inlet  flow  rate  (Anode) 
02  inlet  flow  rate 
N2  inlet  flow  rate 
H20  inlet  flow  rate  (Cathode) 
Cell  current  density 


8.314  [J  mor1  K-1] 

353  [K] 

35 

0.005  [m3] 

0.010  [m3] 

7.034  x  10-4  [mol  s-1] 

7.036  X  10-4  [mol  s"1] 

232  x  10-4  [m2] 

96,485  [A  s  mol-1] 

32  [kPa] 

0.2095 

0.7808 

0.9999 

4.21  x  10-6  [m2  mol  A-1  s"1] 
1.07  x  10-5  [m2  mol  A-1  s-1] 
2 

2.5 

0.8 

0.9 

0.0611  [m3s“1] 

0.0019  [m3  s"1] 

4.5403  [m3  s-1] 

0.1503  [m3  s  1  ] 

0.0019  [m3s-1[ 

100  [A  m~2] 


It  was  found  that  the  rank  of  the  controllability  matrix  is  equal  to 
5,  equal  to  the  order  of  the  system,  n  =  5,  and  hence  this  system  is 
controllable.  In  other  words,  control  inputs  exist  that  can  transfer  the 
state  of  the  system  from  any  location  in  the  state  space  to  any  desired 
location  in  the  state  space  [28],  The  system  observability  is  tested  by 
examining  the  rank  of  the  observability  matrix  defined  by  [28] 


'  C  ‘ 
CA 

O  =  CA2 
CA3 
CA4 


(37) 


It  was  found  that  the  rank  of  the  observability  matrix  is  equal  to 
5  (equal  to  the  order  of  the  system),  and  this  system  is  also 
observable  [28], 


4.2.  Sliding  mode  controller  design 

The  system  in  Equation  (33)  can  be  divided  into  two  subsystems 
as  follow 
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(38) 


<5*,(t)l  [-0.02767  0.00732  0.00732  1  [<5X3(t)"| 

5X2  (t)  =  0.02733  -0.00767  0.02733  <5X4(t) 

dx{t)\  L  0.00005  0.00005  -0.03508 J  |_<5X5(t)J 

[  2.0935  x  10“6  1 
+  2.9127  x  10“3  5Uc(t) 

[  -6.3572  x  10  e  J 
[-9.6493  x  10  61 
+  0  5/(t) 

L  8.1233  x  10“5  J 
=  A28x(t)  +  B25u(t)  +  G28,(t) 

(39) 

The  objective  is  to  keep  the  pressure  differences  of  hydrogen 
and  oxygen  in  a  certain  range  to  protect  the  membrane  damage.  We 
have  designed  the  controller  to  keep  both  pressures  around  3  atm. 

For  the  first  subsystem  of  Equation  (38),  there  exists  a  non¬ 
singular  similarity  transformation  T\  [20],  which  yields 

r$i(t)l  =  [-1.6346  xlO-6  -9.1207  l[qi(t)l 
[q2(t)\  L  8.7110  x  10-7  -8.8370  x  10  3  J  |_q2(t)J 

(40) 

with  d  —  -2.8436  x  10-2.  The  disturbance  can  be  nullified  using 
Drazenovic’s  invariance  condition  since  rank([Bi  Gi[)  —  rank([Bi]) 
[22,29], 

The  pair  (A1mAll2)  is  controllable  since  the  original  system  is 
controllable  [20],  Hence,  we  can  find  a  state  feedback  gain  matrix 
K\  such  that  /\] M  -  K^u  is  asymptotically  stable. 

On  the  sliding  surface  [20],  the  system  trajectory  in  the  (qi  (t),  q2 
(t))  coordinates  is  expressed  as 

[K,  !][’[«]  -  [-0.5482  1){£®]«« 

or 

st(0  —  [0.5482 

#  0.5482x1(t)  +7.9622  x  104x2(t)  =  0 


(41) 


(42) 


in  the  original  coordinates. 

Starting  with  Sr(t)  =  0,  we  design  the  sliding  mode  control  law 
for  the  sliding  surface  (42) 


*i(t>  =  0  =  G,[g]  =  GiAt  [^([)]  +CiB1(u1(t)  +  df(t)) 

(43) 

From  Equation  (43)  and  using  the  result  of  Equation  (42),  control 
Ui  (t)  is  obtained  as 


The  Pressure  of  Hydrogen 


Fig.  2.  Pressure  of  hydrogen  of  the  linearized  system. 


-  -(GiB.r'Cia,  [*’$]  -(C,B,r'{T, +«i) 


where 

Tl  ^jGjBilldW  (45) 

is  required  to  overcome  the  disturbance  /  (t)  and  <ri  >  0  provides 
that 

s1(t)s1(t)<0  (46) 

is  satisfied. 

The  second  subsystem  (39)  does  not  satisfy  Drazenovic’s 
invariance  condition.  Instead  of  using  Utkin  and  Young’s  method 
[20],  we  define  the  sliding  surface  for  the  second  subsystem  as 


s2(t)  =  y2(t)  —y2ref(t) 


(47) 


The  Pressure  of  Oxygen 


Fig.  3.  Pressure  of  oxygen  of  the  linearized  system. 
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Current 


Fig.  4.  Current  changes. 


similar  with  the  method  of  Talj  etal.  [7],  Starting  with  s2(t)  =  0,  we 
design  the  sliding  mode  control  law  for  the  sliding  surface  defined 
in  Equation  (47) 


s2(t)  =  0  =  y2(t)  =  x3(t) 

=  <133*3  (t)  +  ££34*4(0  +  <234*5(0  +  b-i2uc{t)  +g3/(0  (48) 

From  Equation  (48),  control  u2(t)  is  obtained  as 


“2(0  =  -;1  (£233*3(0  +  £234*4(0  +  £234*5(0  +  a2 

D32 

+  |«3|/max|S2(0l  (49) 

where  o2  >  0  is  chosen  to  satisfy  the  condition 


s2s2  <  0  (50) 

Note  that  we  can  use  the  same  procedure  to  find  the  sliding 
surface  and  the  controller  for  the  subsystem  1. 

The  simulation  results  are  presented  in  Figs.  2  and  3.  The 
disturbance  current  waveform  is  presented  in  Fig.  4. 

In  Fig.  2,  the  sliding  mode  control  law  is  from  Utkin  and 
Young’s  method  [20]  and  in  Fig.  3,  the  control  law  is  found  by 
using  s  =  y2  -  y2 ref.  It  can  be  noticed  from  Figs.  2  and  3  that  the 
pressures  are  kept  around  3  atm  after  about  1  s  despite  the  current 
density  changes  during  entire  time  internal  of  interest,  from  0  to 
5  s. 

5.  Conclusion 

We  have  applied  the  sliding  mode  strategy  for  the  linearized 
model  of  the  well-known  nonlinear  model  of  the  proton  exchange 
membrane  fuel  cell  obtained  using  the  MATLAB  Symbolic  Tool  Box. 
For  this  well  defined  model,  which  has  uniqueness  of  steady  state 
variables,  asymptotic  stability,  controllability,  and  observability,  the 
sliding  mode  technique  copes  very  well  with  the  cell  current 


changes  /(t),  and  keeps  very  precisely  the  pressures  of  hydrogen 
and  oxygen  at  the  desired  (required)  values. 


Appendix  A 

vn  =  Cjl2  -  2Ah2C1/(Yh2  -  1  )kaUa  +  A^2Yh2(Yh2  -  2 )k2au2a 
<P12  =  Cjl2  -  Au2ClIYH2kaUa  +  Aj^  (Yh2  -  l)^U2 
<P21  =  VMY h2  -  1  )kaUa  -  Aj^Y2,  -  3Yh2  +2 )k2ti2 
<P22  =  Cjl2  -  Ah2C]  /Yh2 kaUa  +  4  (Yh2  -  l)k2u2 

<P3i  =  I2Cf(  -  C1+c2  +  (C,  +  C2)YN2)  +  Aair/2Ci  (3C: 

-  (C:  +  2C2)YN2)  ( -  4C2Y02  -  2(C,  +  C2)YN2Y02)kcUc 
+  2Xlirl(2Ct  -  Q  YN2  -  4C[  Yo2  +  2(Ca  +  C2 )Y22 

+  (C:  +  2  C2)YN2Y02)k2u2  +  4Aair(YN2  +  Yq2 

-  2)  Y02k2u2Y02 

<P32  =  I3C2(-  C,  +  C2  +  (Ci  +  C2)YN2)  +  Aair/2C! ( 

-  2(3Ci  -  C2)  +  (3C,  +  4C2)YN2  +  2(Ci  +  C2)YQ2)kcuc 
+  2X2akI(2C2YN2  +  (Ci  +  2C2)Y02  -  3Q)k2u2  +  4A^.(YN2 
+  Y02-1)^u^ 

<P4i  =  f2Ci  YN2  ((Ci  -  C2)  -  (Ci  +  C2)YN2) 

+Aair^N2  (  -  3Ci  +  (Cl  +  2C2)YN2)  (2(C,  +  C2)Y02)kcUc 
+2  A2ir(  -  2  +  YN2  +  Yo2)YN2k2u2 

V42  =  A2 Ci  ((C,  -  C2)  -  (Cl  +  C2)  YN2)  +  Aair/YN2  ( -  4Ci 
+  (Ci  +  2C2)YN2)  (2  (Ci  +  C2)Y02)kcuc 
+  2Aair(  —  1  +  Yn2  +  Yo2)^cuc 

V51  =  13C2Yn.2  (Cl  -  C2  -  (Cl  +  C2)YN2)  +  Aair/2C 1  (Bi  +B2  +  B3 
+  B4)kcuc  -  2 Aair/(B5  +  B6  +  B7)k2u2  -  4/.air(2  +  YNj 
+  Y02){\+YN2+Y02)I4u3c 

Bi  =  2 (Ci  -  C2)  -  7CiYN2;  B2  =  -2(Ci  -  C2)Y02 

B3  =  4(Ci  +  C2)YN2  Yq2  ;  B4  =  (3  Ci  +  4C2)Y22 

B5  =  3Ci  +  2 (Ci  -  C2)YN2 ;  B6  =  -(5C,  +  2C2)  Yq2  +  2C2Y22 

B?  =  2(Ci  +C2)Y22  +2(Ci  +2C2)Yn2Y02 

<P52  =  I3C2(-  Cl  +  C2  +  (Cl  +  C2)YN2)  +  Aair/2Ci  ( 

-  2(3Ci  -  C2)  +  (3C,  +  4C2)YN2  +  2(Cr  +  C2)YQ2)kcuc 
+  2A2irJ(2C2YN2  +  (Ci  +  2C2)Yq2  -  3Ca )k^u2  +  4X3ak(Yfi2 
+  Y02-\)k3cu3c 
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Appendix  C 


ail  =  v  (  MhkaUa  ■  Cl!)  — 


ai2  =  y-  {W^Ua  -  Cl 7)  — 


d21  =  - _<PaPvS  2  +  /—  X2-  2)  k“Ua 

VA  \(Xi  +  X2  -  (PaPvs)  (Xl  +  X2)  / 

_Rl(r  ' 

MVl+*2)2j 


b2 1 


+  X2-V«  ^1+^2/ 

b32  =  ^airfy02  -= - _3  _  V 

V  2  x3  +  x4  +  x5;  c 


"  Vc  a,r  V.X3  +  X4  +  X5  -  <paPvs  X3+X4+X5J 


a2  2 


KTX  f  -<PgPvs 

H"  \(Xl  +X2-rpaPvs)2 


033 


RT 

Vc 


Crj\  X4+X5 

2  J  (x3+x4  +  x5)2 


\(x3  +  X4  +  x5  -  <-/)cP„sj2  (X3  +  X4  +  X5j^ 

x3  +x4 
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